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ABSTRACT 

We present observations of Taurus-Auriga Class I/II protostars obtained with 
the Spitzer InfraRed Spectrograph. Detailed spectral fits to the 6 and 15.2 mi- 
cron ice features are made, using publicly-available laboratory data, to constrain 
the molecular composition, abundances, and levels of thermal processing along 
the lines of sight. We provide an inventory of the molecular environments ob- 
served, which have an average composition dominated by water ice with ~12% 
CO2 (abundance relative to H2O), >2-9% CH3OH, -14% NH3, -3% CH4, -2% 
H2CO, -0.6% HCOOH, and -0.5% SO2. We find CO2/H2O ratios nearly equiva- 
lent to those observed in cold clouds and lines of sight toward the galactic center. 
The unidentified 6.8 micron profiles vary from source to source, and it is shown 
to be likely that even combinations of the most common candidates (NH^ and 
CH3OH) are inadequate to explain the feature fully. We discuss correlations 
among SED spectral indices, abundance ratios, and thermally-processed ice frac- 
tions and their implications for CO2 formation and evolution. Comparison of our 
spectral fits with cold molecular cloud sight-lines indicates abundant prestellar 
ice environments made even richer by the radiative effects of protostars. Our 
results add additional constraints and a finer level of detail to current full-scale 
models of protostellar and protoplanetary systems. 
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INTRODUCTION 



The chemical and thermal evolution of protostellar environments affects the spectral 
appearance of ices observed in the sequential stages of pre- main-sequence evolution (e.g. 
Willner et al.l Il982l : iGerakines et al.l Il999l : Ivan Dishoeckl |2006| ) . The presence of the ice ab- 
sorption features is often easily recognized, with the most common ones in the mid-infrared 
wavelength range found at 6.0, 6.8, and 15.2 fim, but the exact feature shapes depend 
strongly on the composition of the ice as determined by the initial condensation and sub- 
sequent thermal processing and UV irradiation. Mixing ices can cause wavelength shifts in 



the resonances of the individual components with in the ice matrix (e.g. lEhrenfreund et al. 
1999l : IPalumbo fc Barattalboool : lOberg et al.ll2007l ). and differences in lattice structure caused 



by either radiation damage or thermal processing are spectroscopically identifiable (e.g. 
Hagen et al.lll98ll : iMoore fc Hudson|[l992l : iMaldoni et al.l[l998l : iKouchi fc Kurodalll99Cl[ ). 



The initial reservoir of star formation material, the diffuse ISM, is extremely ice-poor 



[e.g. iGillett et al.lll975l : IWhittet et al.lll997l). It is not dense enough to shield its dust grains, 
upon which ices condense (e.g. iTielens fc Allamandolalll987l ). from incident destructive pho- 
tons. Only when cold molecular clouds of high density {uh ^ 10^ cm~^) are formed can 
interstellar ices survive and be observed. In the Taurus molecular cloud, for example, a well- 
known seat of starbirth, ices are observed towards background stars even in quiescent regions 



unaff ected by the newly-formed stars (e.g. lBergin et al.ll2005l : iKnez et al.ll2005l : IWhittet et al. 



20071 ). These studies reveal absorptions attributable to H2O, CO, and CO2, as well as some 
less-constrained absorptions (e.g. the still unidentified band at 6.8 /xm) due to presumably 
more complex species. The ice bands observed towards these field stars are smooth, indicat- 
ing amorphous and unprocessed compounds, and low upper l imits can be placed on some 



species created by UV radiation (e.g. "XCN" ; ICibb et al 



20041 ). These observations support 



the theory that molecular cloud ice chemistry is dominated by efficient cold grain surface 
reactions. The prese nce of larger complex molecules (such as HCOOH, which can be formed 
at low temperatures; lKeandl200ll ) demonstrates the potential complexity of these reactions 
and indicates that protostars are born into a pre-enriched chemical environment. 

Surveys of deeply embedded young stellar objects (YSOs) reveal a limited range of 
thermal processing, indicative of differ ences in YSO environment and temperature (e.g. 



Gerakines et al.l Il999 : 
in the classification of 



Gibb et al.l 120041 ). Studies of less-extinguished protostars (Class I 



Ladal 119871 ) show evidence of enhanced complex molecule formation 
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(H2CO, XCN, OCS, PAHs [polycyclic aromatic hydroc arbons]), ice crystallization, further 



thermal annealing, and ice matrix seKregation (see e.g. lAlexander et al 



2OOII : lEhrenfreund et al.lll998l : IWatson et al.ll2004l : iGibb et all 12004 : iBooeert et al.ll2004l ) 



2003 



Keane et al. 



Finally, the least embedded Class II objects have seemingly accreted or dissipated most 
of their cold natal envelopes and have heated the outer layers of their dust disks above 
the subhmation temperature of astrophysical ices. These objects tend t o show silicate dust 



emission, PAH ern issions, and weak or nonexistant ice signatures (e.g. iForrest et al.l 12004 



Furlan et al.ll2006l ). though some Class II edge-on systems have ice abso rptions arising from 



the cold disk midplane (e.g. iPontoppidan et al.ll2005l : iFurlan et al.ll2008l ) 



Differences in chemical and physical struct ure exist between YSOs with masses either 
above or below a few Mq (Ivan Dishoeckl l2003l ). Some of these dissimilarities may be at- 
tributed to observational biases (low-mass YSOs are generally closer and can be better 
resolved), but others are likely to be dependent on the intrinsic nature of the object (e.g. the 
role of shocks, core temperature, and infall timescal es). For example, low- mass YSO studies 



claim higher abundances of CO2, relative to H^^O (jBoogert et al.ll2004j : IPontoppidan et al. 



20081 ) and amorphous sihcate dust (IWatson et al.ll2004j ). than surveys of more massive YSOs. 
This paper focuses on the relatively underexplored categories of low-mass Class I proto- 
stars and Class I/II transition objects (henceforth both groups will be collectively referred 
to as Class I/II objects) with ice ab sorptions by analyzin g mid-infrared spectroscopy of 16 
Taurus- Auriga YSOs (<l-3 Myr old: Iwhite fc Gh^boOll ) , observed wit h the InfraRed Spec - 
trograph (IRS: lHouck et al.ll2004j ) on board the Spitzer Space Telescope ( IWerner et al.ll2004j ). 
Our sample overlaps slightly with that of the extensive recent work by the c2d team (e.g. 



uur sample overlaps siigntiy witn tnat 01 tne extensive recent worK Dy tne cza team [e.g. 
Boogert et aPboOsI : IPontoppidan et allboosl : loberg et al.lboosi ). and we indicate differences 



in our approach and findings throughout this paper. This study helps complete our view 
of the ice evolutionary sequence from cold molecular clouds through warmer environments 
with protoplanetary disks. These low-mass, possibly early solar analogs are essential pieces 
to the puzzle of solar system formation. 

Section [2] contains a short description of the sample, observations, and data reduction. 
In Section [3] we describe how the ice features are identified, isolated, and analyzed, and in 
Section H] we look at trends seen in the feature optical depths. In Section O we derive the 
composition of the carriers of the primary ice absorptions and evaluate candidates for some 
of the less constrained features. In Section [6] we discuss composition and spatial distribution 
of the ices and place our observations in the context of pre-main-sequence evolution. 
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OBSERVATIONS 



2.1. Target Sample Description 



The Taurus- Auriga dark cloud is a nearby (d ~ 140 pc; iBertout et al.lll999l ). low density 



1-10 pc , iGomez et al.lll993l ) star forming region containing predominantly young, 



low-mass protostars. It provides an opportunity to study closely low-mass star formation 
unaffected by the shocks and strong UV radiation associated with more massive protostars. 
Based on IRAS 25 and 60 nm photometry, an exte nsive sample of Class I/II objects was 
selected from the list by iKenyon fc HartmannI (119951 ) and observed with Spitzer's IRS (Sec- 
tion [2]2]). From this sa mple we chose for analysis the subset of 12 Class I and 4 Class II 
objects (as classified in iKenyon fc HartmannI Il995l ) with distinct ice features in regions of 
unsaturated, well-calibrated data. The IRAS spectral energy distribution (SED) points used 
in the original sele ction are insufficient to determine system inclination, but modeling by 
Fur Ian et al.l (120081 ) indicates that our smaller sample spans a wide range of disk inclination 
angles, from nearly edge-on to face-on systems. 

For comparison, we include a secondary sample comprising observations of Taurus and 
Serpens cloud field stars — Elias 16, Elias 13, and CK 2 — as well as of the Class protostars 
CepE-MM and L1448 IRS 2 (AOR IDs 3623680, 3620352). We also include data from the 
line of sight towards the W3(0H) star formation region, which probes part of the foreground 
W3 molecular cloud (AOR ID 3566336). 



2.2. Observational Procedure and Data Reduction 

Our sample was observed with Spitzer's IRS as part of the IRS Disks project (PI: 
J. Houck & D. Watson; PID 2) during guaranteed time awarded to the IRS instrument team 
(2004 February-March). All of our targets were observed using the full 5-37 /im wavelength 
coverage, combining either the Short Low and Long Low (SL, LL; 5.2-14 /im, 14-38 /im; 
A/AA ~ 90) modules on the IRS, or the SL with the Short High and Long High (SH, LH; 
10-19 yum, 19-37 /xm; A/AA ~ 600) modules, depending on the expected fiux levels. 

Most (12) of the objects were observed using a small spectral map (3 x 2) rather than a 
single pointing; this strategy avoids the need for pointing peak-ups and allows, in principle, 
for a better reconstruction of the spectrum in cases of small mispointings. The three-step 
sequence was carried out in the dispersion direction, with the steps separated by three 
quarters (for the SL module) or one half (for the other modules) of the slit width. The 
two-step sequence corresponds to the normal IRS nod distance of one-third the slit length. 
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The other 4 objects were observed in the IRS staring mode, using only the two nods along 
the slit. Observational details of all our Class I /II targets are listed in Table [TJ 



We analyzed data from pipeline version S13, using the smart package (iHigdon et al. 



Furlan et al. 


(2006) and 


, . . 

Furlan et al. 


(2008) 



The uncertainty in each wavelength bin is taken to be half of the nod difference in flux 
level, except where this value is smaller than 1% of the flux. In those cases, the uncertainties 
are taken to be 1% of the flux in order to prevent underestimation. We estimate an absolute 
spectrophotometric accuracy of 5-10%. 



For the field stars, we use the previou sly published data for Elias 16 (IBergin et al.ll2005l ) 
and Elias 13 and CK 2 (IKnez et al.ll2005l ). without re- reducing the data. 



2.3. IRS Spectra Description 

In Figure [H we show the IRS spectra of the 16 Class I/II stars in our sample. The 
YSO mid-infrared spectra are dominated by broad solid state resonances, particularly those 
at 6.0, 6.8, 9.7, and 15.2 fim. Additional, but less pronounced, features are visible in some 
spectra at 7.7 /im and as a broad band centered at 18 fim. Silicates are responsible for the 
broad 9.7 and 18 fim features, which can occur in either absorption (e.g. L1551 IRS 5), 
emission (e.g. IRAS 04108+2803A), or both. 

The 6.0 fim feature is primarily due to the 0-H bending mode of water ice, while the 6.8 
/im feature remains largely unidentified, with CH3OH and NH^ (among others) as suggested 
carriers. The 7.7 fim absorption is generally attributed to solid CH4, and the 15.2 fim band 
is due to the 0-C-O bending mode of solid CO2. A complete inventory of absorption features 
identified in our spectra is contained in Table [2J 

Spectra for W3(0H) and the Class objects (Figure [2]) contain many of the same 
features. 
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3. METHODS 

3.1. Isolating the IR Features and Determining Optical Depths 

3.1.1. YSO Spectra 

The optical depth, r, and shape of the individual absorption features under consideration 
are extracted from the astronomical spectra using a mathematical continuum fit. We make 
a cubic spline fit to the ~5-20 fim continuum of each spectrum using the following relatively 
feature-free regions: bluewards of the 6 /im H2O feature; between the unidentified 6.8 /xm 
feature and the onset of the 9.7 /^m silicate band; between the red wing of the silicate band 
(~14.5 /im) and the 15.2 /im CO2 feature; and redwards of the CO2 feature — while ensuring 
that none of the feature wings are included in the continuum. We calculate the optical depth 
in the spectral features using Fx = F\ cont^^^^- 

Because the 9.7 /im silicate feature is highly complex, and detailed modeling of the 
silicate mineralogy is beyond the scope of this paper, accurate extraction of ice features 
blended with the 9.7 /im silicate complex (e.g. those of CH3OH, NH3, and H2O) is im- 
practical. However, the 11-13 /im H2O libration band and the 18 /im absorption feature 
in particularly silicate-rich spectra are much broader and shallower than the 15.2 /im ice 
absorption they overlap, so we can treat these shapes as pseudo-continua in order to isolate 
the CO2 ice absorption. 

Additional steps are needed to extract small features in the 7-8 /im range. The wide 
variation in 6.8 /im and silicate absorption profile wings requires that each spectrum be 
treated individually. For objects not affected by the silicate feature bluewards of 8 /im, a 
straight-line continuum (~7.1-8 /im) is imposed to isolate any features. For objects affected 
by both the 6.8 and the 9.7 /im features, two straight-line continua are imposed: one from 
~7. 1-7.5 /im and one from ~7.5-8 /im (Figure [3]). 



3.1.2. Laboratory Spectra 



We use laboratory data from the Leiden University ice analog datab asej^ (jCerakines et al. 



I995I . Il996l : lEhrenfreund et al.lll996l . Il999l : Ivan Broekhuizen et al.ll2006l ) to identify the com- 



position and apparent temperatures of the ices seen in our targets. In order to provide a 
direct comparison between the feature optical depths measured in the laboratory and those 



http: / /www.strw.leidenimiv.nl/'^lab/databases/databases.htm 
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observed in the astronomical spectra, we removed a straight-line continuum from the labora- 
tory data. For the wavelength ranges of interest, optical depths calculated using a straight- 
line continuum differ, on average, by less than half of a percent from those calculated using a 
more complicated cubic spline one. The wavelengths of the baseline ends are generally equal 
to the continuum points used in the astronomical spectra's spline fitting, with an occasional 
slight shift to ensure feature wings are excluded from the continuum. In Figure H] we show 
the range of abundances contained in the Leiden databases' non-irradiated H20:CH30H:C02 
mixtures (T=10-185K, all circles); note that this plot does not represent the fraction of the 
database comprising other molecules (e.g. CH4 and H2OH) or UV-irradiated mixtures. 

The temperatures given in this paper, except where otherwise noted, refer to labora- 
tory conditions. Interstellar temperature barriers for crystallization, sublimation, and other 
physical processes are lower than their laboratory counterp arts, due to the low er pressure 
and longer interaction timescales in interstellar space (e.g. iBoogert et al.ll2000l ). Assump- 
tions about particle shapes have not been applied to the laboratory spectra used in this 
broad analysis, and such corre ctions would in any c ase be inappropriate for the annealed 
inhomogeneous mixtures used ( iGerakines et al.lll999l ). 



3.2. Spectral Fitting 



To confirm the validity of using laboratory spectra to model circumstellar material, we 
estimate the total extinc tion Ay along the line of sight towards our YSOs using Ay/rsi ~ 17 
( iRieke fc Lebofskylll985l ). We find Ay ~ 3-34 for our total sample — and Ay ~ 6-34 for those 
objects without clear signs of dust emission contaminating the absorption feature — compare d 
with the average extinction of Ay ~ 1-5 toward Taurus association stars (IMyers et al.lll983l ). 
In addition, recent work has sho wn that in dense regions, this relation with a given rsi 
substantially underestimates Ay ( IWhittet et al.l Il988l : IChiar et al.l 120071 ). so the total cir- 
cumstellar extinction is likely to be even higher. 

Thus, we assume that the foreground absorption towards the Taurus region is rela- 
tively low and that the spectral ice features result from material in t he circumstellar en- 



vironment. In contrast with a least-components type of analysis (e.g. iBoogert et al.l 12008 



Pontoppidan et al.ll2008l ). we determine the ice composition and column densities of the cir- 
cumstellar ice by fitting linear combinations of laboratory data to the observed IRS spectra. 
With the laboratory data described in Section I3.1.2[ detailed spectral fitting is performed 
on the 6.0 and 15.2 fim features using a standard x^-minimization technique. The equation 
solved is 
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where the scale factors of each component are the free parameters in the fit. To compare 
goodness-of-fit among various combinations, we use xt^ the reduced parameter. This 
is defined as "~ '"^); where n is the number of spectral data points and m is the 

number of free parameters in the fit. Traditional statistics require ~ 1 for an acceptable 
fit, but our values are often significantly larger or smaller. Artificially small xl values can 
be caused by uncertainty overestimation or by poor signal-to-noise. Large values may be 
a result of uncertainty underestimation, unknown uncertainties in the laboratory spectra, 
or the existence of components not included in the fits. Therefore xl comparisons may be 
made to determine the best laboratory fit to a single spectrum; however, comparisons of 
goodness-of-fit among different YSO spectra are not necessarily valid. Quantities derived 
using a fitted spectral feature (e.g. processed fractions, column densities) are calculated 
using an average of all fits for that feature with a xt within 20% of the minimum. 



3.3. Column Densities 



[e.g. 



To calculate column d ensities along the lines of sight to the YSOs, we use the relation 

E^aAA 



Alexander et al.ll2003[ ) 



N{cm 



A2 A ' 



(2) 



where AA is the wavelength bin size, Apeafe is the wavelength of the feature's peak optical 
depth, and A (cm molec"^) is the band strength as determined in the laboratory. 



3.4. SED Spectral Index 



To provide better resolution in evolu tionary status than the sini ple iLadal (119871 ) classi- 



fication, we use the spectral index a (e.g. iKenyon fc Hartmannlll995l ): 



dlogXFx 

a = ~ 

dlogX 



log - log AiFa, 
log A2 - log Ai 



(3) 



where we choose Ai and A2 to lie within well-calibrated and nearly-featureless parts of the 
spectrum but as much separated as possible (i.e. Ai = 7 fim and A2 = 30 /im; Table [3]). 
The slope of this section of the spectrum depends upon a variety of factors, including disk 
radius, total envelope density, outflow cavity opening angle, and inclination angle. With the 
exception of inclination angle, all of these factors depend upon the total dust column density 
in such a way that more dust produces a steeper SED slope, and models dem onstrate that 



inclination angle effects are generally less significant than these other factors (IFurlan et al. 
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20081 ). Thus small values for the spectral index a indicate a relatively blue or exposed stellar 



object, but with increasing a the spectrum becomes redder, indicating the YSO is more 
heavily embedded in dust. 



4. OPTICAL DEPTH TRENDS 

In Figure O we show correlations in the optical depths of the three strongest absorptions 
attributed to ices (at 6.0, 6.8, and 15.2 /im) in our YSO sample. Clearly, the 6.0 and 6.8 
yLim features are very closely correlated, with linear Pearson correlation coefficient r = 0.96 
(Figure E^). The peak ratio is Tg.s/Tg.o ~ 0.96±0.03, and the straight-line fit passes within 
0.01 of the origin. Since the 6.0 /zm band is primarily attributed to the bending mode of 
H2O, this correlation indicates that the carrier(s) of the 6.8 /im band forms and evolves in 
a similar, polar environment. 

There is a looser correlation between the 6.0 /zm and the 15.2 yum optical depths in our 
YSOs (r = 0.7, Figure [5)d), resulting in a wider ratio range r6.o/Ti5.2 ~ 0.43-1.65. When the 
same ratios are calculated for the laboratory spectra, though the spread is even wider (0.01- 
3.6), only a small fraction fall within the range defined by the YSOs (Figure HI filled circles); 
the 15.2 /im CO2 feature is generally much stronger than the 6.0 /xm water ice feature, inde- 
pendent of temperature. Mixtures dominated by H2O provide most of the closest matches. 
In addition, there are two lab mixtures of nearly 1:1:1 H20:CH30H:C02 abundances (center 
of Figure H]) which have r6.o/Ti5.2 ratios comparable to the YSO ones but only after they have 
been heated to high temperatures (T > 125 K), nearing the ice sublimation points. The 
Leiden database has a scarcity of water-rich, C02-poor (^20%) H20:CH30H:C02 mixtures 
in the temperature range 10-180 K, though based on this comparison, these compositions 
appear to be astronomically relevant. Because of the undersampling of compositions in this 
range, we are not able to constrain the ice composition well from the T6.0/T15.2 ratio alone. 

Despite being closely correlated with rg.o, the 6.8 /im peak optical depth is even less 
correlated with the 15.2 /im depths (r = 0.56, Figure Et). This indicates that much of the 
inherent scatter in the Tg.s/Te.o ratio is due to variations in the 6.8 /xm feature among different 
lines of sight. 



5. SPECTRAL FEATURE ANALYSES 



In order to set more accurate constraints on the ice properties than feature-strength 
analysis provides, we perform spectral fitting and detailed analysis of the most prominent 
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features. 



5.1. The 6.0 fim Feature 

Each YSO spectrum is fitted over the wavelength range 5.4-6.3 fim using the method 
described in Section [X^ and the component set listed in Table HI An example result is shown 
in Figure O and the total sample residuals, after subtraction of the 6 fim components, are 
shown in FigureO Due to overlap with the unidentified 6.8 fim feature, the range beyond 6.3 
fim was excluded in order to avoid imposing artificial structure on any fit residuals. These fits 
focus on identifying the constituents of the 6.0 fim feature — H2O, H2CO, HCOOH, and NH3, 
as described in the follo wing subsections. This suite of species is similar to that employed 



by lBoogert et al.l (120081 ) to explain their extracted component set Cl-5. 



Either noise or saturation of the absorption bands in the background star/ Class sample 
prevents detailed spectral fitting of the 6 /xm features, so we scaled a pure H2O (10 K) 
spectrum to match the peak optical depth at 6.0 /xm in order to calculate A^(H20) (Table [8]). 
The shape of the 6 fim W3 spectrum is too irregular to be reasonably fit with a solid H2O 
lab spectrum in this manner. 



5.1.1. H2O 

Water is the least volatile and most abundant of the ices typically found along protostel- 
lar lines of sight. The mid-infrared spectrum offers three strong features of water ice: narrow 
ones at 3.1 and 6.0 fim, and a broad (A ~ 11-20 fim) feature peaking at 13 fim. To turn 
these profiles into water ice column density is somewhat problematic in each case. The 3.1 
nm ice feature is g enerally regarded as the cleanest measure of H2O ice column density (e.g. 



Tielens et al.l 1 19841 ). However, inconvenient gaps in the spectrum transmitted by Earth's at- 
mosphere make it difficult to determine the continuum on the short wavelength side of this 
feature in ground-based observations, which would be the source of most of the comparisons 
we could make to our data. In addition, rather significant extinction and scattering near 
3 fim are inevitably associated with protostars and can affect the profile of the water ice 
feature. Futhermore, this may lead to the bulk of the 3 fim background continuum taking 
a significantly different path through the protostellar envelope than that which gives rise to 
the longer-wavelength features. The latter effect would be rather more serious in systems 
viewed closer to edge-on, such as IRAS 04169+2702 and DG Tau B. 



Extinction at 6 fim is similar to that at 13 fim, and both are significantly smaller 
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than that at 3 nm (ex tinction in magnitudes smaller by factors of ~ 2-3; iMathid 11990 
Indebetouw et al.l l2005l ). Thus the optical paths sampled at 6 and 13 /im, and the ice 
column densities derived, should be similar and may be more characteristic of the line of 
sight to the embedded protostar than is the case at 3 /im. In observed spectra, the 6 fim 
feature is the easiest of the three to analyze, as there are fewer problems in determination of 
the continuum. However, the depth and profile of the 6 /im feature are possi bly influenced 



by ab sorption from species other than water ice, as shown, for example, by iBoogert et al. 



(120081 ). The 13 /xm water ice libration feature is broad and overlaps many other absorptions, 
including the loosely-constrained 9.7 /im silicate feature, so its optical depth and feature 
shape are hard to measure precisely. 

Taking account of all of these difficulties, we choose to measure water ice column density 
with the 6 /im feature, and check against the possibility of substantial "contamination" 
from other species using the 13 /im feature. To calculate the H2O column density, we use 
the 6.0 fim feature fit results along with Equation 2 (summed over A=5.5-7 /^m), with a 
band strength of y4ii/2o=l-2xl0^''cm molec"^ for the T = 50 K water ice components and 
74/^20=1-3x10^ ''cm molec"^ for the T = 120 K one. These values pertain to pure H3O ice, 



though they are also suitable for apolar dilutant concentrations of <30% (iGerakines et al. 



I995I ). We note that the derived column densities iV(H20) (Table 5) may in some cases 
represent upper limits, as the overall composit i on of the 6 /im feature is still a source of 
discussion (e.g. Gibb et al. 2004; Boogert et al. 20081 ). However, as we will see below, the 
consistency of these results with those from 13 /im features in our spectra indicates that 
the contribution of unaccounted-for species at 6 /im is rather small compared to that of 
water ice. In general, our values of water ice column density are somewhat larger than those 
derived for the same objects from the 3.1 /im feature (e.g. Boogert et al. 2008), and tend 
to depart more the closer the object is to an edge-on view, so the differences may in part be 
due to the expected differences between the light paths viewed at the two wavelengths. 

To look for signs of additional ice components with features near 6 /im that are not 
resolved spectrally from the water-ice feature, we have measured the 11-20 /im solid H2O 
absorption in the nine targets for which the 10 /im silicate f eature is most deeply in ab- 
sorption. To do so, we followed a procedure similar to that by lBoogert et al.l (120081 ). From 
a low-order polynomial fit to the continuum, constrained to follow the spectrum at 5.3-5.6 
and 30-35 /im, we determined the total optical depth in the 8-30 /im range. To this we 
"fit" th e profile of interstella r silicate absorption, exemplified by the galactic center source 
GCS 3 (IKemper et al.l 120041 ) and constrained to match the observed spectrum at 9.9-10.1 
fim and 23-30 /im. Since this interstellar silicate profile may not be fully applicable to cir- 
cumstellar silicates, this step introduces additional uncertainty. Subtraction of the resulting 
profile from the optical depth spectrum reveals the water ice libration absorption, peaking 
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at ~13 yum, and interrupted at 15.2 /im by the CO2 ice feature. We measured from such 
results the peak optical depth in this water ice feature around 13 /im. By small variation of 
the wavelengths at which the continuum and silicate profile were constrained to follow the 
spectrum, we estimated the uncertainty in the peak optical depth. Because the shape of this 
very broad water feature is more sensitive to both the continuum and silicate profile fits than 
is the peak optical depth, the corresponding water ice column density is considerably more 
uncertain. In Figure [TJ the values for peak optical depth are plotted against the column 
density A^(H20) derived from the 6 /im feature as described above. Evident in Figure [7] is a 
good correlation between the 13 and 6 fim water signatures, with Pearson's r = 0.74. That 
the correlation is tight indicates that, if an additional ice component is significant at 6 /im, 
its abundance relative to water ice must be essentially constant. 



5.1.2. H2CO and HCOOH 



Many of the objects in our sample have additional absorptions centered around ~5.8 
yum that can be well-fit by a combination of H2CO (formaldehyde) an d HCOOH (formic 



acid) ices. Solid H2CO is observed in many protostellar systems (e.g. iKeane et al.l 12001 
Gibb et al.ll2004l) an d plays an important role in the production of complex organic molecules 
(ISchutte et al.lll993l ). H2CO is difficult to detect independently in YSO Spitzer spectra; its 
3.47 fim C-H stretch feature is beyond the wavelength range of the IRS, and its weaker 6.68 
/im C=0 stretch feature is blended with the 6.8 /im absorption, itself not well understood. 
The H2CO column densities (Table [5]) are derived using the observed integrated optical 
depth of the 5.8 /im feature (a s determined from the spectral fits) and a band strength of 
/lj:^^PO=9.6xlO-^^cm molec"^ JSchutte et aDll993h . 



Solid an d gas-phase HCOOH has been detected in both low- and high-m a ss star-forming 



regions fe.e. ICazaux et al. 


2004; 


Keane et al. 


2001 


)• ^ 



Schutte et al.l Il999l : iBottinelli et al.l l2004l : iRemijan et al. 



absorptions at 5.8 and 7.2 fim, and several in the 8-11 /im range (iBisschop et al.ll2007l ). The 
latter ones are blended with the 9.7 /im silicate complex, but the 7.2 /im band, though intrin- 
sically ~4 times weaker than the 5.8 /im resonance, is detectable and indeed is tentatively 
identified in 6 of our 16 objects (Section 15.31) . We have derived HCOOH column densities 
(TableE]) by integrating the observed optical depth of the 5.8 /im feature (as determ ined from 



the s pectral fits) and using a band strength of Ahcooh=Q-7>^^0 ^^cm molec ^ ( jMarechal 
1987h . 
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5.1.3. NH^ 



Solid NH3 (ammonia) has mid-infrared features at 6.15 /im and 9.3 fim (e.g. lSandford &: AUamandola 



19931 ). blended with the 6 /im water ice and 9.7 /xm silicate features, respectivel y. The effect 
on th e 6 fim H2O feature is negligible until the NH3 concentration reaches ~9% (IKeane et al. 



2OOII ). At higher concentrations, an excess appears around 6.15 fim on the red side of the 
6 /im band in laboratory spectra, which fits well the broad 6 fim feature observed in many 
of o ur YSO spe c tra. A nalyses of the 9 /im inversion mode resonance by iLacy et al.l (119981 ) 
and iGibb et al.l (120001 ) indicate that circumprotostellar NH3 exists primarily in an H2O- 
dominated matrix. Therefore, an H20:NH3 (1:0.2) mixture, rather than pure NH3, is used 
in the spectral fits. Because no satisfactory band strength measurements for NH3 in such a 
mixture exist, the upper limit NH3 column densities (Table [5]) are estimated by using the 
band strength of pure water ice for the H20:NH3 mixture and extrapolating the column 
density of NH3 based on its fractional abundance in the final fit. 



5.1.4. PAHs 



After contributions due to H2O, H2CO, HCOOH, and NH3 are subtracted, small resid- 
uals in the 6.1-6.4 fim range, varying in strength, shape, and peak position, are observed in 
11 of the sample's 16 YSOs (Figure [8]). PAH absorption has be en suggested as a ca ndidate 
for the 6.2-6.3 nm feature observed in many massiv e protostars (ISchutte et al.lll996l ). Using 
the average interstellar PAH spectrum compiled by lHony et al. J200lh as a reference, we ex- 
amined the peak position and widths of the YSO residuals. Of the objects within our sample 
with excess absorptions peaking between 6.2 and 6.3 fim, only three (IRAS 04361+2540, DG 
Tau B, and HL Tau) have an excess with a width similar to that seen in the PAH spectrum, 
and none of these three match the PAH feature in shape and peak position. Though the 
PAH relative band strengths often depend upon molecular environment and ionization level, 
the broad 7.5-9 /im complex and the 6.2 /xm feature are co nsistently observed at comparable 
strength in astronomical environments (IHony et al.ll200ll ). yet we find no correlation in the 
YSO spectra between the 6.2 and 7.7 /xm optical depths — the ratios are scattered between 
and 7.5. Moreover, the pror ninent but narrow PAH peak at 11. 3 /im (consistently ~1.4x as 
strong as the 6.2 /im peak; IChan et al.ll2000l : iHony et al.ll200ll ) and the lesser ones at 12.7, 



13.5, and 14.3 /im are not observed in the YSO data. Thus in agreement with iKeane et al. 



(120011 ). we find it unlikely that PAHs contribute significantly around 6.2-6.3 /im. 
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5.2. The 6.8 yum Feature 



Though the 6.8 /xm interstellar absorption band was discovered nearly 30 years ago 
JPuetter et al.lll979h . its identification remains ambiguous and inconclusive. Species ranging 
from hydrocarbons to alcohols t o carbon dust hav e been proposed to contribute to the 
feature in a variety of sources (see lKeane et al.ll200ll ). Here we focus on CH3OH and NH^ — 
the c andidate carriers most commo nly accepted as explaining at least part of the absorption 



le.g- 



Boogert fc Ehrenfreundll2004l ). 



5.2.1. cm OH 



The position of the 6.8 /xm feature coincides with the C-H deformation mode of alcohols, 
and of these, methanol is the most likely to be found in interstellar abundance. Theoret- 
ical models of ice processes predict the creation of solid methanol through the sequential 
hydrogenation of CO ice which has been accreted from the gas phase onto dust grains 
(ITielens fc Charnleyl 119971 ) . The presence of methanol ice in interstellar and protostellar 
environments is observationa lly supported by d etection of its own resonances — e .g. the C- H 



stretc hing mode at 3.54 /zm (IGrim et al.lll99ll ) and the C-0 stretch at 9.8 /im (IGibb et al. 
20001) — as well as by sign atures of its interactions with other molecules (especially CO2; e.g. 
Ehrenfreund et al.lll998l ). 



Our sample spectra do not contain the NIR wavelengths of CHaOH's multiple stretch- 
ing modes, and the 9.8 /im resonance is blended with the silicate absorption. However, 
constraints on the CH3OH abundance can be derived using the 15.2 /xm CO2 band, as the 
interaction of CH3OH with CO2 changes the feature shape near 15.4 /xm (Sections 15.41 and 
16. ip . The relative strength of the 15.4 /xm shoulder implies that at most 50% of the 6.8 /im 
feature is due to CH3OH in a CO2 matrix, although CH3OH unmixed with CO2 could be 
present to provide additional absorption at 6.8 fim. Since CH3OH is thought to form in a 
hydrogenated environment, the assignment of at least part of the 6.8 fim absorption to this 
molecule is then supported by the close correlation of the 6.0 and 6.8 fim optical depths 
(Section m). 



5.2.2. NHt 



The NHJ ion was proposed as the carrier of the 6.8 yum band bv iGrim et al.l ( 



and since then, this proposal has achieved varying levels of success (e.g ISchutte et al. 



989), 



1996 



Keane et al.ll200ll : ISchutte fc Khannall2003l ). The ion can be produced in the laboratory by 
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the acid-base reaction HNCO + NH3 OCN-+ NHj flKeanelllQQTh or in UV-irradiated 
H20:CO:NH3 mixtures (as in Figure [8]). 

NHj has expected counter-ion features at 4.62 fim (OCN~) and 6.3 and 7.41 fim 
(HCOO~, produced by photolysis). The 4.62 fim feature does not fall within the IRS wave- 
length coverage, and the 7.41 fim signature is too weak to be conclusively identified in our 
data. However, some of our sample objects do have a small excess absorption near 6.3 fim 
(Section I5.1.4p . which is unlikely to be attributable to PAHs. It may instead be due to 
HCOO~, indicating photolysis and subsequent warming of a polar ice containing NH3. We 
use the ratios of the laboratory 6.3 fim (HCOO~) and 6.8 /xm (NH4 ) peaks (re.g/Te.s ~ 1.7 
after the H2O feature is removed), along with the equivalent widths of the YSO and labora- 
tory features, to estimate the NH4 contribution to the 6.8 fim feature. Using this method, 
we find that, within our sample, 6-58% of the 6.8 fim feature can be attributed to NH^ (HL 
Tau is an outlier whose particularly large 6.3 /im excess indicates <91% of its 6.8 fim feature 
may be due to NH4 using this method). 



5.2.3. Fits and Summary 

Figure [S] shows a comparison of the 6.8 /xm features, after subtraction of the 6 /xm 
contributors, with a variety of 6.8 yum candidate laboratory spectra, all rebinned to the 
resolution of the observations. The bottom two lab spectra (Lab B and C) have CH3OH/H2O 
abundance ratios of 100 and 0.1, and the gray vertical lines indicate the peak positions of 
these mixtures. The relative amount of methanol present affects the peak substructure. 
While in some cases the YSO peak positions and relative strengths closely resemble the 
laboratory data (e.g. IC 2087 IR and IRAS 04239+2436), in others there is no correlation 
at all (e.g. IRAS 04169+2702). The abundance of CH3OH (relative to H2O) has been found 
to vary by almost an o r der of magnitude among different protostellar lines of sight (e.g 



Pontoppidan et al.l l2003l : iGibb et al.l 120041 ) , which may partly explain the observed variety 



of peak substructures. 

The top laboratory spectrum in Figure [8] (Lab A) is a H20:CO:NH3 deposit, irradiated 
to create NH4 and subsequently warmed from 10 to 165 K. The irradiated spectrum at 
10 K peaks 0.1 /im bluewards of our sample's bluest feature; in order to align with any of 
the observed profiles, the irradiated mixture must be heated to near the H2 O sublimation 



temperature, as the sublimation appears to cause the redward peak shift (ISchutte et al. 



19961 ). Moreover, the width of the NH4 feature is too narrow (~0.25 jiw) to account for 
the entire 6.8 /xm feature widths observed towards each YSO (average ~0.4 /xm). Even if 
multiple NH^ laboratory spectra at various temperatures are combined in order to reproduce 
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the observed feature width, the resultant peak position remains too blue to match the YSO 
data. 

Noting that many of the YSO 6.8 fim features share characteristics with the profiles 
of both CH3OH (triple peaks, 7 /im shoulder) and NH4 (asymmetric peak near 6.75 /im), 
we attempted fitting them with a combination of the laboratory components. While no 
remarkably good fits were achieved, there emerged a general trend that equal contributions 
from NH^ and from CH3OH in a water-dominated matrix, with varying amounts of nearly- 
pure CH3OH intermixed, came closest to reproducing many of the YSO feature width-to- 
depth ratios (i.e. mixtures of irradiated NH3, CH3OH/H2O~0.1, and CH3OH/H2O~100 in 
the approximate ratios 1:1:0 to 1:1:1). However, the overall 6.8 /xm astronomical feature 
shapes could not be be well-matched in terms of either peak alignment or subpeak relative 
strengths. 

Based on this work, we conclude that we cannot rule out NH4 and CH3OH as contrib- 
utors each of at most ~50% to the 6.8 /xm band, but we emphasize that even the 50/50 
combinations do not reproduce the overall feature s hapes. Studies of t h is feature performed 



with higher-resolution ISO-SWS observations (e.g. iKeane et al.l I2OOII : ISchutte et al.l Il996l ) 
find a similar variety of structure and peak positions, confirming that it may be difficult or 
physically impossible to adequately explain the feature using the same set of components for 
all sight-lines. 



5.3. The 7-8 /im Range 

Because of the uncertainties due to the low resolution and low signal-to-noise of the 7-8 
fim optical depth spectrum, as well as the inherent narrowness of the 7-8 /im features, we 
did not attempt to fit the shape of the spectral features with laboratory profiles but rather 
focused on detection of the 7.25 fim HCOOH, ~7.6 fim SO2, and 7.7 /xm CH4 features. 



5.3.1. HCOOH 



The 7.25 /xm H COOH band is intrinsically ~4 times weaker than the 5.85 fim one 



(iBisschop et al.l 120071 ). Based on the amount of HCOOH in the fitted spectra in Section [5?H 
the expected 7.25 /im optical depths in the YSO spectra are <0.02. Indeed, we do observe 
small features at 7.25 /xm in 6 of our 16 objects. It should be noted, however, that the strength 
of these do not correlate with the HCOOH column densities derived from the spectral fits and 
the stronger 5.85 /xm feature (r = 0.06); there are also other species with weak signatures 
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very close to this wavelength (e.g. HCOO^ and HCONH2; ISchutte et al.lll999l ). As the 
optical depth profiles here are neither strong nor resolved enough to conclusively distinguish 
among these species' signatures, these HCOOH detections must remain tentative. 



5.3.2. CH4 and SO2 



In 13 of our 16 objects, we observe a feature near 7.7 /im, where the (z/4) C-H deformation 
mode of CH4 is expected. Laboratory studies have shown that this resonance's position and 
width are strongly dependent on the molecular environment, becoming blue-shifted (as far 
as 7.64 /xm) and b roadened when mix ed with a variety of other ices, including H2O, CO, 
CH3OH, and CO2 (IBoogert et al.lll997l ). In addition, there is the possibility of blending with 
the nearby SO2 feature, whose 7.58 /im S-0 stretching band has been shown to vary its peak 
position from 7.45-7 .62 /xm, depending on the temperature and the molecular environment 
JSoogert et al.lll997h . 



In order to disentangle absorption due to SO2 and CH4, we fit each YSO's 7.45-8 yum 
wavelength range with two Gaussian curves, allowing the peak of one (representing SO2) to 
shift between 7.45-7.62 /xm and constraining the other peak (CH4) to 7.62 /xm or larger. A 
Gaussian approximation was chosen because both the SO2 and CH4 profiles are too narrow, 
compared with the IRS resolution, to evaluate any substructure. In many cases, a single 
Gaussain was sufficient, indicating a single, distinct feature of either SO2 or CH4, depending 
on peak position. Of our 16 YSOs, 7 show only a 7.7 /im feature with little or no asymmetrical 
excess, 3 only have a distinct feature peaking within the SO2 range, 4 show both distinct 
features, and 2 show a broad absorption with additional peaks within the SO2 range (i.e. 
fit by overlapp i ng G aussians). Using the range of laboratory feature FWHMs measured by 
Boogert et al.l (119971 ). we find that in 4 of the 9 objects with either distinct features or a 
Gaussian fit in the 7.45-7.62 /im range, the feature width is too narrow (AA <0.06 /im) to 
be attributed to solid SO2. These may be due to random noise creating an apparent peak or 
to CH4 gas fin es found between ^7 .4-8 /im. SO2 also has gas lines in this region, but these 
are very weak (IBoogert et al.lll997l ). For the 5 r emaining features, w e use Equation 13.31 with 
a band strength A5O2=3.4xl0~^''cm molec"^ (IKhanna et al.l 119881 ) to determine the SO2 
column density; for the other objects, an upper limit to A^(S02) is determined (Table [5]). 

Of the 13 objects with a distinct 7.7 /im CH4 feat ure or Gaussian fit, 7 have FWHMs 
comparable to those observed in laboratory mixtures (IBoogert et al.l 119971 ). particularly in 
those mixtures where CH4 and either H2O, CII3OH, or NH3 exist in roughly equal abundance. 
The remaining 6 objects have wider features, which may be explained by variations in tem- 
perature or ice mantle composition along the line of sight. To determine the column density 
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of CH4 towards each YSO, we us e Equation 13.31 with a band strength 74(7ff^=7.3xlO~^^cm 
molec"^ (pure CH4; discussion in iBoogert et al.l 119971 ) and integrate the optical depth be- 
tween 7.6-7.9 fim. Upper hmit estimates for the 3 non-detections are provided by the noise 
in the spectrum (Tabled]). 



5.4. The 15.2 fim Feature 



Double-peaked structure in the 15.2 fim feature, often observed in systems with an 
exposed source of heating and radiation (such as Class I/II objects), is associated with the 
dual peaks of the degenerate z/2 bending mode in pure CO2 ice. This structure becomes less 
pronounced or disappears altogether when the apo lar CO2 molecules are embedded in an 
H20-rich polar matrix (e.g. lEhrenfreund et al.lll999l ). Thus the presence of a double-peaked 
feature indicates that thermal processing has occured, segregating the CO2 from the polar 
matrix. The broad shoulder on the red side of the absorption, centered around 15.4 /im, 
has been attributed to CO2 interactions with other molecules, particularly CH3OH, and the 
strength and shape of the shoulder depend on both temperature and relative abundances (e.g 
Dartois et al.lll9991). Theoretica l studies suggest that so lid CO2 is formed by oxidation of CO 



Tielens fc Hagenlll982l : iRuffle fc Herbst 



20011) ■ so the co-ex istence and interaction 



ice (e.g. 

with solid CH3OH (also formed from CO; iTielens &: Charnley 1997) is not unexpected. 

In the fits, made for the high- resolution spectral^ between 14.7 and 16 /xm, we include 
cold (10-30 K) deposits of H20:C02, pure CO2, and an irradiated deposit of C0:H20, 
as well as warmer (>100 K) mixtures of H20:CH30H:C02 (Table [6]). The cold mixtures 
containing H2O are used to estimate the contribution from nonprocessed material along the 
line of sight. The warmer H20:CH30H:C02 spectra represent ices which have been warmed 
and have undergone the chemical and physical changes, including segregation of the CO2, 
resulting from that processing. The pure CO2 component, though measured at 10 K, does 
not greatly differ in shape until sublimation at ~80 K and is here used as a measure of 
segregated, highly-processed ice. 



Our fitting components differ from the recent work by iPontoppidan et al.l (120081 ). which 
is focused almost exclusively on detailed methanol-free CO:C02:H20 mixtures. Certainly 
these combinations can produce good fits to the CO2 ice feature. However, in the follow- 
ing we will adhere to the conventional approach involving a warm CH3OH ice component, 
since we find that in most cases this produces better fits to our spectra than methanol-free 
mixtures, and it takes into account the strong evidence of solid methanol in many protostel- 



^The low-resolution spectra are insufficient to resolve any structure which might be present. 
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lar environments (e.g. iPontoppidan et al.ll2003l : iGibb et al.ll2004j ). It is worth noting that, 
since good fits are obtained by both classes of models, neither method should be thought to 
provide a unique solution to the ice-composition analysis. 

Using a single component fit, we find that a nearly equal-parts mixture of H20:CH30H:C02 
at temperatures of >100 K provides the lowest-x^ fit to the observed YSO 1 5.2 /xm absorption 



bands, in accordance with p revious results for more massive protostars (e.g. lEhrenfreund et al. 



19991 : iGerakines et al.lll999l ). Where present, the dual peak positions align, and the approxi- 
mate strength of the 15.4 /zm shoulder is replicated. The fit improves, however, by an average 
of ~50% of Xu^ when we allow multiple components (Figure [9]). The combination of these 
additional components reproduces the relative peak strengths, the inter-peak dip depth, and 
the shape of the 15.4 fim shoulder without requiring a specific single temperature. 

Our multicomponent fits indicate that for our YSO sample, a smaller amount of CH3OH 
(~14%) at lower temperatures is sufficient to account for the observed 15.4 /im shoulder, 
and equal parts H2O, CO2, and CH.sOH are refiected at high temperatures (higher than 
those identified by e.g lEhrenfreund et al.lll999l ). when the CO2 has at least partly migrated 
out of the matrix. Laboratory mixtures at these higher temperatures that also contain the 
smaller fractions of CH3OH have an asymmetrical dual peak that is not observed in the 
YSO spectra; mixtures with H2O, CH3OH, and CO2 in a 1:1:1 ratio at lower temperatures 
show a stronger shoulder at 15.4 /xm (relative to the main peak) as well as a nearly-single 
peak at 15.2 /im. While the cold-surface formation processes and expected abundanc es of 



CH3OH are not tightly-constrained (e.g. Tielens fc Whittet 19961 Keane et al. 2001 . and 



references therein), the generally-higher abundances observed towards high- mass protostars 
are evidence that warmer environments may be conducive to additional CH3OH production. 
This is consistent with the stronger 15.4 /im shoulders observed towards massive p rotostars 
(compared with low- mass ones; e.g. lEhrenfreund et al.lll998l : iBoogert et al.ll2004l ). and the 
equal-parts H20:CH30H:C02 mixtures required to fit them. Most of our YSO fits are further 
improved by including the very cold (~30 K) C0:H20 layer, representing absorption along 
the line of sight due to cold material, affected by the region's radiation field but too distant 
or too shielded from the central source to have undergone thermal heating. 

Methanol-b earing ice mixtures have distinctive absorptions elsewhere in the mid-infrared 
spectrum, notably at 6.8 and 9.8 /im. The former feature is discussed above (Section l5.2p . 
The latter unfortunately lies in the very bottom of the silicate absorption feature. In the 
mixtures that best fit the 15.2 /im ice feature, the 9.8 /im feature is either very weak or not 
particularly sharp. Thus it is difficult to separate possible 9.8 /xm methanol ice absorption 
from small silicate-feature profile variation. In general, a combination of absorption from the 
ice mixture that best fits the 15.2 /xm feature, and from interstellar silicate grains (toward 
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GCS 3; iKemper et al.ll2004l ). can be made to fit the spectrum, as illustrated in Figure [TOl 
However, this result depends sensitively on the fit of the continuum and the silicate absorp- 
tion profile, which as noted above is not very well-defined for these types of environments. 
We found no meaningful, quantitative constraints on a putative CH3OH ice component from 
the 9.8 fim feature. 

The 15.2 fim features in the molecular cloud sight-lines have pr e viously been fitte d with 
predominantly cold, polar H20:C02 mixtures (IBergin et al.l l2005l : iKnez et al.l l2005l ). We 
find that the xl of fits can be matched or reduced by including components containing 
methanol — specifically, the cold counterparts of the H20:CH30H:C02 components used in 
our Class I/II fits described above fFi gurefTTl). Th e perc entag e of polar ice fou nd (~75- 
80%) is comparable to that identified by lBergin et al.l (120051 ) and lKnez et al.l (120051 ). and the 
abundance ratios A^(C02)/iV(H20) from those studies are preserved. 

We have calculated the total CO2 column density (Table [5]) observed along each embedded- 
object and molecular cloud line of sight using Equation [2l The exact band strengths of some 
of the laboratory mixtures used in the 15.2 /im fits are undetermined. The fitting results 
indicate that the vast majority of the CO2 ice resides alongside H2O ice, so we use a band 
strength of 1.5xl0~ ^^cm molec~^, which is appropriate for the 15.2 fim absorption from a 
H20:C02 1.6:1 ice (iGerakines et al.lll995l ). The observed optical depth spectra were inte- 
grated between 14.7-16 fim. We use the relative contributions of each fitting component 
to estimate the column density of solid CO2 in each stage of thermal processing (Table [7]): 
radiation-only (10-30 K, "cold"), some thermal heating ("warm"), and CO2 segregation. We 
have also estimated the column density of solid CH3OH embedded in CO2 using the the 
CH3OH abundance and relative contribution of each laboratory fit component (Table [5]). 



DISCUSSION 



6.1. Ice Composition: Abundances With Respect to H2O 



6.1.1. CO2 Column Density 



The column densities observed in our sample of Class I/II YSOs yield A^(C02)/A^(H20) 
~ 0.12(±0.04) (Table Figure [T^ . which is consistent with the values pr eviously derived 
for many different lines-of-sight, averaging iV(C02)/-/V(H20) ^ .17(=b0.03') (IGerakines et al. 



1999; Nummehn et al. 2001 



Gibb et al. 



2004; Knez et al. 2005; Alexander et al. 



2003h. but 



significantly lower th an th e value of N {CO 2) IN {^20 ) ~ 0.3, with large scatter, reported by 
Boogert et al.l (120041 ) and iPontoppidan et al.l (120081 ). Our study is unique in the sense that 
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we use the band strength Aco2 of an H20:C02 mixture. If we instead use the band strength 
for pure CO2, we find A^(C02)/A^(H20) ~ 0.16(±0.04), nearly equivalent to the results of 
the majority of previous studies, which include high- and low-mass YSOs, background field 
stars sampling molecular clouds, and galactic center sight-lines sampling diffuse interstellar 
material and molecular clouds. 

The dashed line in Figure H] identifies the locus of compositions with CO2/H2O = 
0.12, which are clearly underrepresented in the heated samples from the Leiden Univer- 
sity database. Good profile fits can only be achieved by combining ice mixtures, and this 
gap in composition-space may explain some of the difficulty experienced in fitting 6 nm fea- 
tures. Further laboratory exploration of the water-rich mixtures around the observed ratio 
with CO2 may help to improve constraints on the column density of H2O ice based on the 6 
fim feature. 



6.1.2. Evolutionary Trends in CO2 



When considering the Class I/II sample in conjunction with the Class objects, we 
find A^(H20) to be well-correlated (A^(C02) less so) with the spectral index a, where more 
embedded sources (higher a) show higher column densities, as is expected (Figure \TE\ top 2 
panels). The ratio A^(C02)/A^(H20) increases as the objects evolve towards lower a- values 
(Figure [T3| bottom panel). 



The sublimation temperature of solid H2O is ~40 K higher than that of pure CO2 ice 
(IGerakines et al.lll995l ). Thus with increased exposure to thermal radiation, A^(C02) / A^(H20) 
is expected to go down, seemingly contradicting our observations. Although CO2 molecules 
embedded in a polar matrix may become "trapped" and able to survive to temperatures 
nearing the sublimation limit of H2O, the thermal destruction of H2O before CO2 is unlikely. 

An explanation of the high column density of CO2 in more exposed environments may 
be found in the formation of additional CO2 in Class I/II objects. Photodissociation of 
H2O leads to the production of rad icals (e.g. OH) which c an react to form CO2 on the 



grain surface (e.g. C0+0H^C02; lAUamandola et al.lll988l ). Some other reaction path- 



ways to form CO2 in molecular clouds and star-forming regions have ac tivation barriers 



which become less obstructive in w armer, more exposed environments (e.g lRuffle &: Herbst 



2OOII : Ivan Dishoeck fc Blakdll998f ). Indeed, in UV-irradiated H20:CH30H:C02 laboratory 
samples, CO2 eventually re aches a constant abundan ce level due to an equilibrium between 
destruction and formation (lEhrenfreund et al.lll999l ). while the abundance of H2O further 
decreases due to photodissociation. 
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The values measured towards the background stars Elias 16, Elias 13, and CK 2 cannot 
be extrapolated from this trend. The A^(C02)/A^(H20) ratios for these lines-of-sight (Ta- 
ble [H]) are not lower than the ratios observed towards the YSOs. This discrepancy may be 
due to the fact that the interstellar lines-of-sight probe a variety of environments and compo- 
sitional structure, while protostellar sight-lines are dominated by the circumstellar material. 
Possibly, the lower temperatures in interstellar regions stifle CO2 formation mechanisms with 
activation barriers, while areas of increased density favor H2O production. 



6.1.3. HCOOH 



^^cm ^ (L1551 IRS 5), the column densities 



Apart from one outlier at iV(H20) ~ 1.1 x 10 
of HCOOH and H2O are strongly related (r = 0.61, Figure [H]). This supports laboratory 
experiments demonstrating that the most likely formation environments of solid HCOOH are 
those that also produ ce H2O ice, through either radiative processing or cold surface reactions 
jBisschop et al.lboOTh . We find A^(HC00H)/A^(H20) ~ 0-2%, with an average of 0.6%, well 



within the <5% abundance rat i o observed in lines- of-sight towards massive and low-mass 
protostars (jBisschop et al.l 120071 : iBoogert et al.ll2008l ). 



6.1.4. CH3OH 

The observed A^(CH30H)/A^(H20) ratios in our sample do not show a clear trend 
(Figure [H]). It is known that the A^(CH30H )/A^(H20) abundance varies greatly from source 



to source (e.g lBoogert &: Ehrenfreundll2004l ). and it may be better to analyze the CH3OH 



abundance with respect to CO2. 

CH3OH can be detected in the spectrum in three places. When it is intermixed with 
CO2, a shoulder at 15.4 fim due to CH3OH can be observed in the 15.2 fim feature; and 
CH3OH also shows 6.8 and 9.8 fim bands (the latter is blended in the 9.7 fim silicate complex). 
In Figure [15] we show the ratio of the optical depth at 15.4 fim (vertical axis) and 6.8 fim 
(horizontal axis), both with respect to the optical depth at 15.2 yum. Filled circles represent 
the astronomical measurements, while the plus signs indicate the ratios measured in the 
laboratory spectra. The laboratory spectra show a correlation between the two features, 
whose exact position depends on the temperature. The 15.4/15.2 /xm r ratios observed in 
the astronomical data are within the range measured in the laboratory, but the corresponding 
6.8 fim optical depths are significantly greater than the laboratory values. This could be 
understood if only a small fraction of the CH3OH along the line of sight is embedded in a CO2 
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matrix or if there are significant contributions to the 6.8 fim feature from other components. 



Using the laboratory fits to the 15.2 fim complex in the high spectral resolution data, we 
estimate the column density of solid CH3OH intermixed with CO2 as 30-56% with respect to 
CO2, corresponding to A^(CH30H)co2/iV(H20) ~ 2-9% (Table [5]). Assuming that <50% of 
the 6.8 /xm feature is due to to CH3OH (Section [5.2.31) . we estimate iV(CH30H)totai/-/V(H20) 
< 5-17%. There are two outlying values, for IRAS 04108+2803A and IRAS 04181+2654B, at 
32 and 28%. Although we do not place high confidence in these as absolute values because 
of the 6.8 /xm uncertainties, they are consistent with the detection of solid CH3OH with 
abundances (relative to H2O) in t his range up to ~25% towards low-mass YSO s (as well as 
the overall scatter in abundance; iPontoppidan et al.ll2003l : iBoogert et al.ll2008l ). Thus our 
data support the presence of a significant C02-poor H20:CH30H ice component. 



6.1.5. CHa 



We observe a range of A^(CH4)/A^(H20) ratios from ~0.3 to 8% in our sample (Fig- 
ure [H]), with an average of 2.6% and a single outlier (IRAS 04154-^2823) at 23% . Our 
measurements largely overlap with the ratios of ~0.5-4% reported in the lit erature (e.g. 



Boogert et al.lll996l : iGibb et al.ll2004l ). though ratios up to 13% are reported by lOberg et al. 



( 120081 ). The higher values in our sample (>4%, dashed line in Figure [TBI) may be partially 
ascribed to addit ional absorptions un resolvable at the SL-module resolution, such as gas- 
phase CH4 lines (iBoogert et al.l 119971 ). or to matrix-dependent variations in band strength. 
Hudgins et al.l (119931 ) found that the band strength for the 7.7 /xm CH4 feature may vary by 
as much as a factor of 2.8 between H20-rich and CO-rich matrices. We use the value for 
pure CH4, which lies approximately midway between the two extremes. 

Laboratory studies of solid CH4 indic ate that it can be formed by a variety of processes — 
including UV photol ysis of CH3OH ices (jAUamandola et al.lll988l ) and grain-surface hydro- 
genation of carbon (iTielens fc Hagenlll982l ) — and that the resulting a bsorption profile ha s 
characteristics unique to its formation pathway and embedding matrix (jBoogert et al. 1997). 
Unfortunately, the low signal-to-noise of the CH4 detections, and the low spectral resolution, 
prevent a detailed analysis of the spectral profiles; however, the broad range of peak positions 
and widths implies that the CH4 embedding matrices vary across the sample. 
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6.2. Segregation of CO2 Ice 



In the colder parts (<20 K) of a star-forming region, well-removed from the hot cores 
and embedded protostars or buried in dense clumps, ap olar CO2 ice is formed embedded in 
a polar H2O matrix (e.g. iBoogert &: Ehrenfreundl 120041 ) . The 15.2 fxm CO2 resonance will 



appear to be broad and single-peaked. As these ices are heated, as for instance happens 
in the circumstellar environments of YSOs, the CO2 can migrate out of the polar matrix 
and aggregate, either as a sur f ace layer or in bulk apo lar clumps within the matrix (e.g. 
Sandford fc AUamandolal Il990l : lEhrenfreund et al.l Il999l ) . The spectral appearance of the 
CO2 resonance changes, and a double-peaked structure, with peaks at 15.1 and 15.25 /xm, 
becomes apparent rather rapidly (Figure [T7I) . The migration process also destabilizes the 
C02:CH30H interactions and reduces the 15.4 /xm shoulder. Due to the apparent spectral 
difference between segregated CO2 ice and CO2 ice embedded in an H20:CH30H matrix, it 
is possible to determine the relative mass contained in both of these phases. 

In Figure [TS|, we show the variation in the CO2 15.2 /xm absorption feature in our 
sample, along with laboratory spectra of CO2 fully embedded in the matrix (top; 10 K) and 
fully segregated CO2 (bottom; 130 K). The dashed lines indicate the positions of the double- 
peaked segregated CO2 and the shoulder at 15.4 /im due to embedded CO2 interacting with 
CH3OH. We show only the Class I/II objects for which high resolution spectra are available, 
as the double-peaked structure is not resolved in the low resolution spectra. Several of 
the spectra clearly show the double-peaked str ucture, while a l l devi a te from the smooth , 
unp rocesed Taurus clou d profiles as published by IWhittet et al.l ( 120071 ) , iBergin et al.l (120051 ) , 



and iKnez et al.l (120051 ). Indeed, from the multi-component fits described in Section 15.4 
we find evidence in all of our low-mass YSO spectra for moderate-temper ature thermal 



processing (<100 K; CO2 warmed but not s egreg ated), in agreeme nt with iBoogert et al. 



( I2OO4J ) but in contrast with iNummelin et al.l (120011 ). However, unlike IBoogert et al. 
we also find evidence for high-temperature processing — 5 of the 9 objects (IRAS 04016-1-2610, 
IRAS 04239+2436, IRAS 04365+2535, DG Tau B, and IC 2087 IR) have 15.2 /im profiles 
entirely dominated by warmed ices in which a significant fraction (~20-40%) of the CO2 has 
segregated from the polar formation matrix. 

From the multi-component fits, we can derive the mass distribution of CO2 ice over 
the phases (cold, warm, segregated). The phase transitions correspond to laboratory tem- 
peratures, which in turn translate to (lower) astrophysical temperatures. The onset of CO2 
segregation occurs at a temperature near 105 K in the la boratory (Figure [T71) . which cor- 
responds to an astrophysical temperature of ~65 K (e.g. IBoogert et al.ll2000l) . The CO2 
evaporation temperature depends on the relative abundance of other species, particularly 
H2O, but the mixtures with the most physically likely abundances (i.e. CO2/H2O < 1) 



- 25 - 



show an average sublimation temperature of ~150 K, or 80 K in astrophysical environments. 
Thus, the segregated column densities of Table [7] represent CO2 ice at circumstellar radii 
where the ice temperature is between 65 and 80 K, and the cold+warm column densities in- 
dicate ice farther from the central source (i.e. where T < 65 K). The inner boundary of the 
segregated CO2 ice zone is defined by the evaporation t e mper ature at ~80 K. The low-mass 



protostellar envelope models produced by iFurlan et al.l (120081 ) for nearly all the stars in our 
sample indicate these approximate temperature ranges in the densest inner regions {R < 
100 AU); however, detailed physical structure models that also include the appropriate ice 
opacities are needed to further constrain the location and processing levels of the ices. 

In Figure [19] we show the correlation of thermal processing levels with the SED spectral 
index (indicative of the line-of-sight dust column). While there is a large amount of scatter 
in these plots, the overall fraction of cold CO2 ice increases (and the warm fraction decreases) 
for more deeply embedded objects, while the segregated fraction does not show a clear trend 
with total dust opacity. 



6.2.1. Segregation of CO2 Ice in Molecular Clouds and W3(0H) 



All three interstellar cloud sight-lines show single-peaked 15.2 /xm profiles with very 
little substructure; CK 2 and Elias 16 are remarkably similar in shape (Figure [TT]) . despite 
indi cations that ice niantles in Serpens may be fu ndamentally different from those in Taurus 
(e.g lEiroa fc HodappI Il989l : IWhittet et al.l 120071 ). These results give support for the cold 
grain origin of at least some of the CH3OH observed towards protostars; the presence of the 
unidentified 6.8 /xm band in background star sight-lines (in particular CK 2 and Elias 16; 
Knez et al.ll2005l ) is consistent with the assignment of CH3OH as a carrier. 



Interestingly, the spectrum of W3(0H), which has the reddest SED of our entire sam- 
ple, shows a clearly double-peaked CO2 feature, indicating pure CO2 along the line of sight, 
which could have been formed by thermal segregation (Figure I2U1) . The segregation fraction 
is ~26%, comparable to the Class I/II average of ~30%. Given the difference in segregation 
with the pure molecular cloud lines of sight, the most likely explanation is that this line 
of sight samples not only the cold foreground W3 cloud but also the warmer, thermally- 
processed CO2 ices around the massive protostars in the W3(0H) region. The strong sim- 
ilarity of the 15.2 fim absorption profile towards W3(0II) to those of the low-mass Taurus 
objects is curious, as the line of sight toward W3(0H) is highly divers e, containing an 
ultra - compact Hll regiori , strong outflows, a nd H2O/OH maser emissions (ITurner fc Welch 
1984 : lAlcolea et al.lll993l : lArgon et al.ll2003l ). Future spectroscopic observations with higher 
angular resolution will resolve the spatial structure and provide useful information on the 
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complexity of the early stages of CO2 heating in embedded objects. 

7. CONCLUSIONS 

We have analyzed the mid-infrared molecular ice features in 16 Class I/II low-mass 
YSOs using Spitzer IRS observations. Using detailed comparisons to laboratory ice analog 
spectra, we find the ice mantle compositions to be dominated by H2O and CO2 ice but with 
significant abundances of additional species such as CH3OH, H2CO, HCOOH, NH3, CH4, 
and SO2 (<10% each). 

Detailed spectral fits to the 6 //m feature reveal that >80% of the band can be attributed 
to H2O, with additional absorptions near 5.85 //m identified as HCOOH and H2CO. Some 
of the observed spectra have a broader 6 //m peak that can be well-fit by H20:NH3 ices, 
providing evidence for the presence of solid NH3. On the long- wavelength wing of the 6 
/xm feature, small absorption remnants are observed, varying in strength, shape, and peak 
position among the sources. We show that these are highly unlikely to be the result of PAH 
absorption, due to the absence of other expected PAH features at 7.5-9 and 11-14 fj,m. 

We have examined the unidentified 6.8 /xm feature and discussed its assignment to 
solid CH3OH and NH^, the two most common candidates. While no high-quality spec- 
tral fits could be achieved using these two components, we find that on average the ob- 
served YSO features could be approximately fit by equal contributions from NH^ and from 
a CH3OH/H2O~0.1 mixture, with varying amounts of nearly-pure CH3OH intermixed. How- 
ever, we find a wide variety of feature shapes in our sample, indicating that the interstellar 
component(s) producing this band vary greatly from object to object, and a single simple 
carrier or set of carriers may be inadequate to explain the feature. 

In the 7-8 /xm range, we are able to isolate and identify features attributed to CH4, SO2, 
and (tentatively) HCOOH ices. The high CH4 ice abundances derived from these, along 
with the variety of feature shapes and peak positions within the sample, suggest variations 
in line-of-sight molecular environment or possible absorption contamination by unresolvable 
gas-phase CH4 lines. The solid SO2 and CH4 band profiles are very sensitive to molecular 
environment and temperature, and additional high-resolution studies of this wavelength 
range will yield further information on molecular ice interactions and abundance structures 
along the line of sight. 

The 15.2 //m solid CO2 bending mode feature is well-reproduced by mult i- component 
laboratory fits comprising ice analogs (H2O, CH3OH, CO, CO2) in various stages of radiative 
and thermal processing. Using the relative contribution of the components to each high- 
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resolution feature profile, we find that all of the YSO ice environments have undergone at 
least some thermal heating, and 5 of the 9 have profiles entirely dominated by processed ices 
in which a significant fraction (~20-40%) of the CO2 has segregated from the polar formation 
matrix. 

The molecular ice abundance ratios, particularly A'"(C02)/A'"(H20), are consistent with 
the values observed in a wide variety of sight-lines, but within the sample, wc observe 
a negative correlation between this ratio and the mid-infrared SED spectral index. This 
evidence describes a more complicated polar CO2 environment than indicated by previous 
studies. Comparison of the SED spectral index with the ice temperatures derived from the 
15.2 /xm fits confirms the presence of cold ices in more deeply embedded objects and warmer, 
annealed ices in more evolved ones. 

Comparison with molecular cloud sight-lines show that our fits are consistent with a 
scenario in which H20-rich grain mantles with CO2 and CH3OH relative abundances of 
~20 and ^10%, respectively, as well as traces of other species, form in the cold molecular 
clouds and are subsequently heated by young protostars. This thermal processing sublimates 
the most volatile species while inducing molecular interactions and migrations that become 
spectroscopically detectable. The wide range of both molecular abundances and thermal 
processing levels observed within our sample indicates that young stars have a strong effect 
on their surrounding environments. 

We are grateful to Edwin Bergin and Claudia Knez for providing us with the field 
star spectra, to Melissa McClure for assistance with data reduction issues, and to Adwin 
Boogcrt for stimulating discussion. Extensive use was made of Leiden University's Sackler 
Laboratory ice analog databases and of the NASA ADS Abstract Service. This work is 
based on observations made with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, Cahfornia Institute of Technology under NASA contract 1407. 
Support for this work was provided by NASA through Contract Number 1257184 issued by 
JPL/Caltech. E.F. acknowledges support from a NASA Postdoctoral Program Fellowship. 
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Table 1. Overview of the Spitzer IRS observations of our sample of Class I/II objects. 
The first column gives the object name, and the second to fourth columns give the 
observations' AOR key, date, and IRS modules/mode used (SL: short-low, SH: short-high, 

LL: long- low, LH: long-high). 



Object 


AOR 


Obs. Date (2004) 


IRS Modules/Mode 


IRAS 04016+2610 


3528960 


2/07 


SL,SH,LH Map 


IRAS 04108+2803A 


3539712 


2/06 


SL,LL Stare 


IRAS 04108+2803B 


3529472 


2/06 


SL,SH,LH Map 


IRAS 04154+2823 


3534336 


3/03 


SL,LL Map 


IRAS 04169+2702 


3534848 


3/04 


SL,LL Map 


IRAS 04181+2654A 


3546112 


2/29 


SL,LL Stare 


IRAS 04181+2654B 


3545856 


2/07 


SL,LL Stare 


IRAS 04239+2436 


3530752 


2/08 


SL,SH,LH Map 


IRAS 04295+2251 


3537408 


2/27 


SL,LL Map 


IRAS 04361+2547 


3533056 


2/29 


SL,SH,LH Map 


IRAS 04365+2535 


3533312 


2/27 


SL,SH,LH Map 


IRAS 04381+2540 


3538944 


2/27 


SL,LL Map 


DG Tau B 


3540992 


2/08 


SL,SH,LH Stare 


HL Tau 


3531776 


3/04 


SL,SH,LH Map 


IC 2087 IR 


3533312 


2/27 


SL,SH,LH Map 


L1551 IRS 5 


3531776 


3/04 


SL,SH,LH Map 
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Table 2. Absorption features identified in our YSO spectra. 



Wavelength (fim) 


Carrier(s) 


Mode 


A (lO^^^cm molec-i) 


Reference 


Ice Features 


5.81 


H2CO 


C=0 stretching 


0.96 


1 


5.85 


HCOOH 


C=0 stretching 


6.7 


2 


6.0 


H2O 


H-O-H bending 


1.2-1.3^ 


3 


6.8 


CH3OH?, NHJ?, organics 


O-H bending, C-H deforming 






7.58 


SO2 


S-O stretching 


3.4 


4 


7.7 


CH4 


C-H deforming 


0.73 


4 


15.2 


CO2 


O-C-O bending 


1.5'' 


3 


Other Features 


9.7 


Sihcates 


Si-O stretching 






14.97 


CO2 


gas-phase 1/2 bend 




5 


18-20 


Sihcates 


O-Si-O bending 







depends on temperature; range represents 10-120K 
for H20:C02 mix, see Section lOl 



References. — (1) ISchutte et ai] ^99^, (2) iMarechall l|l987h . (3) iGerakines et al.l ^99^, (4) iBoogert et al] {Tgol), (5) 
Ivan Dishoeck et ahl lll996h 
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Table 3. Characteristics of our sample of Taurus Class I/II objects. The first column gives 
the object name, the second and third columns give the 6.0 /im and 15.2 fim peak optical 
depths, and the fourth column gives the SED spectral index a (7-30 /im; Equation [3]). 



Object Te.o^m Ti5.2^m SED a 

IRAS 04016+2610 0.38(±0.01) 0.73(±0.01) 0.43(±0.02) 

IRAS 04108+2803A 0.19(±0.01) 0.19(±0.01) 0.29(±0.10) 

IRAS 04108+2803B 0.26(±0.01) 0.43(±0.01) 0.21(±0.07) 

IRAS 04154+2823 0.08(±0.01) 0.16(±0.01) -0.35(±0.04) 

IRAS 04169+2702 0.25(±0.01) 0.31(±0.01) 0.98(±0.02) 

IRAS 04181+2654A 0.2(±0.01) 0.3(±0.01) 0.08(±0.02) 

IRAS 04181+2654B 0.32(±0.03) 0.44(±0.01) 0.51(±0.03) 

IRAS 04239+2436 0.25(±0.01) 0.44(±0.01) 0.41(±0.05) 

IRAS 04295+2251 0.2(±0.02) 0.18(±0.01) 0.47(±0.02) 

IRAS 04361+2547 0.45(±0.01) 0.27(±0.03) 2.22(±0.02) 

IRAS 04365+2535 0.33(±0.01) 0.63(±0.01) 0.54(±0.04) 

IRAS 04381+2540 0.44(±0.01) 0.62(±0.01) 0.61(±0.03) 

DG Tau B 0.34(±0.04) 0.3(±0.01) 0.38(±0.07) 

HL Tau 0.12(±0.01) 0.12(±0.01) 0.34(±0.02) 

IC 2087 IR 0.14(±0.01) 0.33(±0.01) -0.89+0.01) 

L1551 IRS 5 0.54(±0.01) 0.58(±0.03) 1.41(±0.02) 
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Table 4. Laboratory spectra used in 6 //m fits. 



Composition 


Temperature (K) 


Radiation^ 


Pure H2O 


10 




Pure H2O 


30 




Pure H2O 


50 




Pure H2O 


120 




Pure H2CO 


10 




Pure HCOOH 


10 




H20:CH4^ 1:0.33 


50 


1 hr 


H20:NH3 1:0.2 


10 


1 hr 



^The equivalent of lO^^cm ^ photons (E^ > 6 
eV) for the noted length of time 

does not have a resonance at G/im., but its 
presence in a water ice matrix affects the shape of 
the feature. 



References. — All spectra from Leiden Univer- 
sity's Sackler Laboratory ice analog databases. 



Table 5. Derived solid-state molecular abundances in our YSO sample. The first column gives the object name, and 
the second through eighth columns give the column densities of H2O, CO2, H2CO, HCOOH, NH3, SO2, CH4, and 

CH3OH. 



Object 


A^(H20)^ 




iV(H2C0) 


iV(HCOOH)'= 


Ar(NH3) 


A^(S02) 


iV(CH4) 


iV(CH30H)^ 




lO^^cm^^ 


% H2O 


% H2O 


% H2O 


% H2O 


% H2O 


% H2O 


% H2O 


TT? A Q H/ini fiU-OKl fl 


c: 01 -1-0 9 


1 7 /I _L 1 ^ 


1 Q 
1 .y 


1 9. 
1.0 


^19 1 
<. IZ. 1 


<.V.l 


0/1-1-0/1 
U.4 It U.4t 


0.0 It 1.4 


TT? A Q n/i 1 ns-LOsn"? A 

ixlAo U41Uo-t-ZoUO-A- 


9 on -1- n "? 


-u 9 /I 


u.o 


U.D 




17-1-11 
i . ( It i . 1 


s n -1- 9 A 

c5.U It Z.D 




IrtAo U4iUo+ZoUoiJ 


A A A 4-0 9 
4.44 it V.Z 


1 1 1 _L 1 C 
1 1 .1 zt 1.0 





U.o 


<y.u 


<u. / 


c; A -U 1 9 
O.D it i.Z 


D.o it i.O 




1 n7 -1- n 9 
1 .U ( It u.z 


1 Q _L 7 Q 

ly .0 =t ( .J 


n 1 


u.u 


^9S A 




9*^ 9 -1- Q 1 
Zo.Z zn y . i 




IRAS 04169+2702 


4.33 ± 0.2 


10.1 ± 1.5 


2.4 


1.0 


<9.0 


<0.5° 


3.9 ± 1.2 




IRAS 04181+2654A 


3.63 ± 0.2 


11.4 ± 1.8 


3.6 


0.2 


<9.0 


<o.r 


3.1 ± 1.4 




IRAS 04181+2654B 


4.18 ± 0.4 


15.6 ± 2.3 


7.1 


0.7 


<25.6 


<0.3 


4.5 ± 2.0 




IRAS 04239+2436 


4.44 ± 0.2 


12.6 ± 1.5 


2.4 


0.0 


<12.5 


0.6 ± 0.2 


<0.5 


3.8 ± 1.5 


IRAS 04295+2251 


2.82 ± 0.3 


9.4 ± 2.3 


1.7 


0.4 


<15.3 


<1.0'= 


2.4 ± 1.8 




IRAS 04361+2547 


7.38 ± 0.3 


4.9 ± 1.5 


0.0 


0.9 


<9.0 


<0.1 


4.3 ± 0.7 


1.7 ± 1.5 


IRAS 04365+2535 


5.31 ± 0.2 


15.2 ± 1.6 


1.8 


0.5 


<13.0 


<0.1 


0.8 ± 0.7 


8.5 ± 1.6 


IRAS 04381+2540 


7.92 ± 0.2 


11.8 ± 0.8 


0.8 


0.5 


<9.2 


0.2 ± 0.1 


0.3 ± 0.3 




DGTauB 


5.65 ± 0.6 


7.5 ± 1.5 


0.0 


0.5 


<17.7 


1.2 ± 0.6 


<0.4 


3.5 ± 1.5 


HLTau 


1.91 ± 0.2 


7.3 ± 2.8 


5.9 


0.7 


<11.0 


1.0 ± 0.5 


<1.2 


4.0 ± 2.8 


IC2087IR 


2.13 ± 0.2 


19.5 ± 4.5 


0.8 


0.1 


<26.7 


<0.2 


1.8 ± 2.2 


8.7 ± 4.5 


L1551IRS5 


10.90 ± 0.2 


8.0 ± 1.3 


0.0 


0.0 


<9.0 


<0.0 


1.7 ± 0.4 


3.2 ± 1.3 



Note. — iV(H2C0) and 7V(IIC00II) have typical absolute uncertainties of --10%. 

'^Uncertainties listed are statistical, and derived from the residual differences between the data and our models. There are potential 
systematic errors as well, as described in Section r5. 1.11 



''Using band strength of H20:C02 mixture; see Section WTH 
'^Based on 5.8 /im feature alone. 

'^iV(CH30H) intermixed with CO2, from high-resolution 15.2 //m feature fits. 

''A feature was detected in the SO2 wavelength range but with an overly-narrow width; value given is an upper limit. 



00 

00 
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Table 6. Laboratory spectra used in 15 /xm fits. 



Composition Temperature (K) Radiation^ Reference 



H20:CH30H:C02 0.9:0.3:1 98 ••• 1 

H20:CH30H:C02 1:1:1 130 •■• 1 

C0:H20 10:1 30 1 hr 2 

H20:C02 10:1 10 •■■ 3 

Pure CO2 10 • • • 3 



^The equivalent of 10^^ cm ^ photons (E^ > 6 eV) for the noted length 
of time 



References. — (1) lEhrenfreund et al.l (|l999l ). (2) 
(3) lEhrenfreund et al.l (|l99d l 



Gerakines et al, 



(|l995l l. 
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Table 7. Percentage and column density N (lO^^cm^^) of CO2 ice present in each stage of 
processing along the line of sight: cold, warm, and segregated (due to rounding, 
percentages may not sum to exactly 100%). 



Cold Warm Segregated 



Object 


% Total 


N 


% Total 


N 


% Total 


N 


IRAS 04016+2610 


3 


0.28 


71 


7.3 


26 


2.6 


IRAS 04108+2803B 


33 


1.6 


30 


1.5 


37 


1.8 


IRAS 04239+2436 


10 


0.57 


87 


4.9 


2 


0.14 


IRAS 04361+2547 


62 


2.2 


12 


0.45 


26 


0.95 


IRAS 04365+2535 


2 


0.17 


61 


4.9 


37 


2.9 


DG Tau B 


1 


0.05 


76 


3.2 


23 


1.0 


HL Tau 


32 


0.45 


32 


0.44 


36 


0.51 


IC 2087 IR 


9 


0.35 


68 


2.8 


23 


0.95 


L1551 IRS 5 


22 


1.9 


54 


4.6 


24 


2.1 
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Table 8. Molecular cloud and Class abundances. 





iV(H20) (IQi^cm-^) 


A^(H20) (10i**cni-2) 


iV(C02)/iV(H20)^ 


A^(C02)/A^(H20) 


Source 


[This work] 


[Lit., (Refs)] 


[This work] 


[Lit., (Refs)] 


CK 2 [Serpens] 


2.9 


3.5 (1) 


0.29 


0.33 (2) 


Elias 13 [Taurus] 


0.95 


1.0 (3) 


0.14 


0.15,0.19 (2,5) 


Elias 16 [Taurus] 


2.2 


2.5 (4) 


0.2 


0.18,0.23 (3,5) 


CepE-MM 


<29.7'' 




>0.09 




L1448 IRS2 


<28.6^ 




>0.1 




W3(0H) + W3 











^N{C02) calculated using the band strength A of pure CO2 in order to facilitate comparison with published 
values. 

''Calculated from a scaled pure H2O spectrum. 
References. 



1 



Eiroa fc HodappI |l989l ). (2) lKnez et all |2005h . (3) lwhittet et al.l |l988h . (4) Ichiar et al 



(|l995l ). (5) iNummelin et al 



( 20011 



r 
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Fig. 1. — Our sample IRS spectra. The lower left object has been marked with the locations 
and primary carriers of the most prominent absorption features. Note that the 9.7 /xm silicate 
feature may also appear partly or wholly in emission. 
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Fig. 2.— IRS spectra of the Class protostars CepE-MM (in Cepheus) and L1448 LRS 2 
(in Perseus), and the hne of sight to W3(0H) through the W3 molecular cloud. The 9.7 /xm 
silicate absorption feature in CepE-MM is saturated. 
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Wavelength (/i.m) 



Fig. 3. — Example of the applied continuum determination using the spectrum of L1551 IRS 
5. The inset shows the additional continuum imposed on the 7-8 /im range. 
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Fig. 4. — Range of non-irradiated H20:CH30H:C02 laboratory compositions for which there 
is infrared spectroscopy in the Leiden University database (all circles). The vertices repre- 
sent the pure molecular ices, and distance from each molecule's vertex indicates decreasing 
relative abundance of that molecule (solid lines at increments of 20%). The filled circles 
represent mixtures whose r6,o/Ti5.2 ratio falls within the range spanned by the ratios from 
the YSO spectra (Section H]), and the dashed line indicates the compositional locus where 
CO2/H2O=0.12 (Section EH). 
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Fig. 5. — Relationships between feature peak optical depths for our low-mass YSOs. Panel 
(a) shows Te.s^m vs re.o^m, panel (b) shows ri5.2,„« vs Te.o^m, and panel (c) shows ri5.2^m vs 
Te.s/xm- The dashed lines indicate the best linear fit. 
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5.5 6.0 6.5 7.0 

Wavelength (/xm) 



Fig. 6.— The 6 //m feature of IRAS 04016+2610, fit by a combination of pure H2O at 10 

and 50 K (dots and short dashes), pure H2CO at 10 K (dot-dash), pure HCOOH at 10 K 
(dash-triple dots), and H20:NH3 10:2 at 10 K (long dashes). The solid line is the sum. The 
6.8 /xm feature is deliberately excluded from the fit. 
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1 peak r 



Fig. 7. — Peak optical depth in the water ice hbration mode (A ~ 13 yum), plotted against 
the water ice column density iV(H20) derived from the 6 /im water ice feature as described 
in Section 15.1. 1[ Points and uncertainties are results for the nine Taurus objects in which 
the 9.7 fim silicate features are strongly in absorption. 
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T — I — I — I — I — I — I — I — I — I — I — rrr 




Wavelength (/xm) 

Fig. 8. — The sample's 6-7 fim profiles (after subtraction of the H2O contributions (dashed) 
and all contributions (sohd) at 6 /^m) normalized to 1 and shifted vertically, sorted by 6.8 
fim peak position (the point near 5.6 /im indicates the approximate optical depth uncer- 
tainty). The 5.6-5.7 /im residuals are almost entirely byproducts of continuum mismatches. 
Also shown for comparison are normalized laboratory spectra — Lab A: H20:CO:NH3 10:2.5:2 
at 165 K (irradiated to create NH+), Lab B: H20:CH30H:C02 9:1:2 at 10 K, Lab C: 
H20:CH30H 1:100 at 10 K. 
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1 — ^ 




14.8 15.0 15.2 15.4 15.6 15.8 16.0 
Wavelength (fim) 



Fig. 9. — 15.2 //m feature of IRAS 04365+2535, shown repeated with a vertical offset. 
The spectrum on the bottom is overlaid with a fit of HaOiCHgOHrCOa 0.9:0.3:1 at 98 K 
(dots), H20:CH30H:C02 1:1:1 at 130 K (dashes), and C0:H20 10:1 at 30 K (dash-dot). 

The solid line is the sum. The heavy dashed fit on the upper spectrum is a fit using only 
H20:CH30H:C02 1:1:1 at 118 K. The excess between 14.95 and 15 //m is due to gas-phase 
CO2 ro-vibrational 1/2 Q-branch lines, blended together at this resolution. 
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Fig. 10. — Model of the absorption around 10 yum in IRAS 04 016+2610. Shown a re the 
observed data (points), the GCS 3 sihcate profile (long dashes; iKemper et al.ll2004l ). pure 
water ice at 10 K (dash-triple dots), and the two methanol-bearing mixtures described in 
Section 15.41 and listed in Table [6] (short dashes and dot-dashes), scaled per this object's 
15.2 /im fit. The remnant peak at 9 /im is likely NH3 ice. The model and observation are 
consistent with one another in the relevant wavelength range. 
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Fig. 11. — Comparison of fits to the Serpens background star CK 2 (top) and Taurus 



backg round star Elias 16 (bottom). Tfie left panels are repeats of the fits by iKnez et al. 



(120051 ) ■ where the dotted line is a laboratory spectrum of H20:C02 10:1 at 10 K, the dashed 
line is HsOiCOs 1:1 at 10 K, and the dash-dot line is COiNsiCOa 10:5:2 at 30 K. The 
mismatch on the red wing is caused by the forced continuum at ~16.5 fim, which ignores the 
extent of some of the ice feature wings as measured in the laboratory. In the right panels, 
the dash-dot line is the same, but the H20:C02 10:1 spectrum has been replaced with 
H20:CH30H:C02 9:2:1 at 30 K (long dashes) and the H20:C02 1:1 with H20:CH30H:C02 
0.9:0.3:1 at 30 K (dash-triple dot). 
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Fig. 12.— Correlation between solid iVCH^ O) and iVCCQ^) 



Our sample comprises low-mass protostars. iGerakines et al. 



protostars and some galactic center sources, iNummelin et al. 



or ma ny different lines of sight. 
(ll99S ) contains mostly massive 



(120011) comprises low- and high- 



mass YSOs as well as some molecular cloud sight-lines, and iBoogert et al.l (120041 ) contains 
two low-mass protostars. The solid line is the ratio between column densities using only 
the current work's sample, and t he da shed line is the ratio using all objects except the 
upper outlier from lBoogert et al.l (120041 ). Note that this work uses a CO2 band strength (at 
15.2yum) ~36% higher than that of the previous studies, which may account for most of the 
difference in N{C02)/N{'R20). 
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Fig. 13. — Relationships between the SED spectral index (a; 7-30 //m) and A'"(C02), 
A^(H20), and A^(C02)/Ar(H20). The filled circles represent our Class I/II data, and the 
open circles represent the Class objects. 
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Fig. 14.— Ratios between N{B.20), and A^(HCOOH) and A^(CH30H), for the Class I/II 
sample (Section [6?II) . The bottom plot refers to A^(CH30H) embedded in a H20:C02 matrix. 
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Fig. 15. — The correlation between 6.8 /xm and 15.4 //m optical depths in the laboratory 
ice spectra (plus signs) and the observed YSO spectra (solid circles). Both have been scaled 
by peak 15.2 /im optical depth. The laboratory data cover temperatures from 10-185K 
and CH3OH/CO2 abundance ratios of 0-1.4, which are roughly indicated by the solid lines: 
[<0.5] / [0.5-1.4] / [^1.4]. The vertical spread depends largely on temperature, with warmer 
mixtures of a constant composition showing weaker 15.4 //.m shoulders. 
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Fig. 16. — Comparison of solid A^(H20) and iV(CH4) for the YSO sample. The dashed 
line indicate s the maximum o f ~4% (CH4 relative to H2O) observed towards massive pro- 
tostars fe.g. iGibb et al.ll2004l ). and the dotted line is the newer 13% maximum reported by 



Oberg et al.l (120081 ) towards low-mass protostars. 
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14.50 15.25 16.00 50 100 150 200 

Wavelength (/xm) Mixture Temperature (K) 



Fig. 17. — The "dip-to-peak" ratio as a measure of the amount of CO2 segregation, (a) — 
Demonstrates how the values were measured for each laboratory ice spectrum, using pure 
CO2 as an example, (b) — The relative dip strength as a function of laboratory ice tempera- 
ture, showing the rapid onset of the process near ~105 K. Much of the scatter in (b) is due 
to the dependence of CO2 segregation temperature on the exact composition of the matrix 
in which it is embedded. The large asterisks in the upper left represent pure CO2 ice. 
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14.8 15.0 15.2 15.4 15.6 15.8 16.0 16.2 
Wavelength (/xm) 

Fig. 18. — The high- resolution 15 /im Class I/II profiles visually ordered by "peakiness," 
normalized to 1, and vertically offset for clarity. The dotted fines at 15.1 and 15.25 //m mark 
the peak positions of pure CO2, and the one at 15.4 iim indicates the CH30H:C02-complex 
shoulder. The top and bottom spectra are laboratory mixtures as indicated. 
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Fig. 19. — Relationship between the SED spectral index (a; 7-30 /xm) and the CO2 processing 
fractions in our Class I/II sample. 
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Fig. 20. — 15.2 fxni CO2 feature probing W3(0H) and the foreground W3 cloud. It 

has been fitted with a combination of H20:CH30H:C02 0.9:0.3:0.1 at 98 K (dotted hne), 
H20:CH30H:C02 1:1:1 at 130 K (dashed hne), C0:H20 at 10:1 at 30 K (dot-dashed hne), 
and H20:C02 10:1 at 10 K (dash-triple dot). The solid line is the sum. 



